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The adiabatic theory of resonance one-phonon Raman scattering in semiconductor nanocrystals 
is revised and extended with perturbative non-adiabatic corrections, given by the Albrecht's B 
term. This theory is confronted with the time-dependent perturbation approach, pointing at their 
differences and similarities, ft is shown that both theories are equivalent in the limit of weak 
electron-phonon coupling and non-degenerate or uncoupled resonant states. Evaluations of the 
A and B terms for the confined LO phonon in CdSe and CdS nanocrystals are reported. These 
evaluations show that the B term can usually be neglected. 



I. INTRODUCTION 

In the interpretation of resonant Raman scattering in 
quantum dots, two theories have been mostly used. The 
first one is based on the Albrecht's theorytl of light scat- 
tering from small molecules, where the excited levels that 
contribute to the Raman polarizability are considered as 
vibron states in the adiabatic approximation. This the- 
ory also considered perturbative non-adiabatic correc- 
tions, but ttwft'n JMftfi have been neglected in quantum 
dot studiesBQuoB'UPij. Moreover, the exciton-phonon 
coupling factor has been identified with the Huang-Rhys 
factor, using this as a fitting parameter for the relative 
intensities of different orders in Raman spectra. How- 
ever, the calculations of the BLuang-Rhys, parameter for 
intrinsic exciton states in PbSo and CdSeESO nanocrys- 
tals have given values that are too small compared with 
those needed to explain the experimental results of mul- 
tiphonon Raman scattering. Extrinsic mechpaisms such 
as donpt-iike excitonE3, surface hole trapsEj or extra 
chargestJt-a have been invoked to resolve the discrepancy 

A different model of Raman scattering in quantum dots 
was conceived from a solid state point of view. In this 
approach, the Raman cross sections are calculated from 
third or hig hc i j Q f4er time-dependent perturbation the- 
ory (TDPT)ll3il3tLZl and the intermediate virtual states 
are considered as tensor products of electronic states, 
lattice vibrations, and photons. A non-perturbative cal- 
culation of multiahonon Raman spectra have been re- 
cently presented!!.! Not having adjustable parameters 
this model has scarcely been used by experimentalists to 
interpret their data. Moreover, up to now it is unclear the 
relation between the TDPT and the Albrecht's theory. 

The purpose of the present article is: (1) to establish 
the relation between Albrecht's theory and TDPT, and 
(2) to investigate the importance of non-adiabatic correc- 
tions within Albrecht's theory. The structure of the pa- 
per is as fallows. First, we give an overview of both theo- 
retical approaches and show their interrelation. Next, we 
explain the calculation of the Albrecht's A and B terms 
for semiconductor nanocrystals. Finally we discuss the 
numerical results obtained for several types of nanocrys- 



tals and present our conclusions. Several mathematical 
steps are given in the appendixes. 



II. THE THEORIES 

A. Time-dependent perturbation theory 

For a one-phoncpjRaman process the differential cross 
section is given byE-3 
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where r\i {tj s ) is the refraction index at the incident (scat- 
tered) light frequency uii {co s ), V is a normalization vol- 
ume of the radiation field, and c is the velocity of light in 
vacuum. The one-phonon transition amplitude Mpj {p) 
can be calculated by time dependent perturbation theory, 
considering the unperturbed Hamiltonian as the sum of 
the Electronic, Lattice, and Radiation operators 

H = He + Hl + Hr, 

while the perturbation is the sum of Electron-Lattice, 
Electron-Radiation, and Lattice-Radiation interactions 
(the last one is negligible in resonance conditions) 

Hint = He-l + He-r + Hl-R. 
The electron-lattice interaction can be expressed as 

He-l 



+ 
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where D n > and b v (6J,) are annihilation (cre- 

ation) operators of electronic and vibrational excita- 
tions, respectively, while Hi_ L {v) and its Hermitian ad- 
joint Hp_ L {v) are operators that act on the electronic 
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systemta. The operators He-r and Hr-l have a similar 
structure to that of H E _ L . The phonon states created by 
b\ are those of the electronic ground state G. The elec- 
tronic excitations fj, ^ G can be considered as electron- 
hole pairs and confined excitons. Hence, the one-phonon 
Raman transition amplitude can be calculated by third 
order perturbation theory as 



E 



(hwi - E m ~ huj p + iT^) 
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where T M are the lifetime broadenings of the elec- 
tronic excitations. The interaction matrix elements 

(Ma |#e_l(p)| Mi), (G /i 2 ) and (mi |-^e-k| ^) 

have been calculated under two approximations of solid 
state theory: the Effective Mass Approximation and a 
Long Wave Continuous Model for the optical phonons 
(a field theory fiy T Ypa,GhA ..Wpijlring expressions can be 
found elsewhereffiO'BBBEI The conditions of va- 
lidity of the above formalism can be summarized as: (1) 
The scattering process is dominated by extended vibra- 
tional states, which are not affected by single-electron ex- 
citations, as is usual in large molecules and solids; (2) Ex- 
cited electronic states are well separated in energy from 
the ground state. 



B. Albrecht's theory 

In this scheme, the wave functions of the molecule (or 
the quantum dot) are considered in the adiabatic Born- 



Oppenheimer approximation 



*«.({r}, {<?}) = e e ({r}, {Q})<M{Q}), 



where e and v are the sets of electronic and vibrational 
quantum numbers, respectively; {r} is the ensemble of 
electron coordinates of the molecule (or the quantum 
dot), and {Q} is the ensemble of normal coordinates of 
the ions. The vibrational wave function $ e „ is factored as 
a product of the wave functions of all the normal modes 



<M{Q}) = <(QiK 2 (Q 2 )..., 



(4) 



v a being phonon occupation numbers. According to the 
dispersion theory, Albrecht obtained the Raman polariz- 
ability tensor for resonance scattering as 



a = A'" + B + h.o.t.+n.r.t. 



where h.o.t. means higher order terms and n.r.t. non- 
resonant terms. 



The A and B terms are given by 
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In the above expression we follow Albrecht's notationS, 
but we use dyadic notation instead of tensor subscripts, 
fkoi instead of hfo for incoming photon energy, and 
f)w a p = E a — Ejj. \gi), \ev), \gj) are the initial, interme- 
diate and final vibrational states (Eq. @)), respectively, 
in the potential energy fields of the electronic states g 
(ground) and e (excited). 9Jt° g and h® e are, respectively, 
dipole and electron-phonon interaction matrix elements 
(see Table |). Also, the order of some matrix element 
indexes is reversed in order to generalize the Albrecht's 
expressions for the case of complex wave functions. Note 
that the order of subscripts in matrix elements in this 
notation is the opposite of that in Dirac's notation. The 
summation on the index (e) can be restricted to the reso- 



nant state e. The s states in (pf) appear from a perturba- 
tive expansion of 6 e ({r}, {Q}) in terms of Q s ({r}, {0}). 

The above formalism is valid if: (1) The e states are 
non-degenerate or uncoupled to other states with the 
same energy, and (2) Excited electronic states are well 
separated in energy from the ground state. 



1. The offset oscillators model 

In this model, the vibrational states ipf (Q a ) and 
<P% (Qa) are assumed to be localized in parabolic 
Vg!e)(Qa) potentials with the same curvature, but the 
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TABLE I: Equivalence of Albrecht and TDPT notations. 



Magnitude 



Albrecht 's 



TDPT 



Electronic 
excited states 
Electronic 
ground state 

Matrix element 

Matrix element 



e, s 



G 



^ ei -m° a!S (^\H-_ R ( ei ,0)\G) 
xoKe <M2 \He-l(p)\vi) 



Raman spectra. The fitted values of axe-near 1, in 
contradiction with microscopic calculations&JO. 



2. The limit of weak electron-phonon coupling 

In the limit of weak electron-phonon coupling the os- 
cillator offset A a should be small. Expanding ([?]) and (||) 
in powers of A a we find 



origin of V e shifted in AQ a — \/2x a0 A a = \f2^/hjuj a A a . 
In many cases, only one vibrational mode is assumed 
and A a is employed as fitting parameter. Moreover, A 2 a 
is identified as the Huang-Rhys factorBEEl. The integrals 
appearing in the A and B terms are given by 



(9l\ev) = J^e-^A^'L^ (A») 



(7) 



and 



(ev\Q a \gl) = x a0 



l + l 



(g,l+l\ev) 



+ x aO\J g {g,l~ Mev) 
l-v + A 2 a/1 . . 

Xa0 ~7^A — ^ 9 
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where L v m are the Laguerre generalized polynomials and 

XaO = y/h/uj a . 

Several authors have used this m °dcl|-to n study the 
electron-phonon coupling in nanocrystalsa&tMtjtJ. All of 
them have considered only the term A and have used 
A^ to fit the overtone/fundamental intensity ratios of 



(ev\gi) = 5 v>i + 
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— x a0 
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i\ I v\ 

' —r ((i + l)8 v i + S v i-2) - \ -7T (i$v,i + 8 v ,i+2) 



Replacing (g) in Eq. (||) we obtain 



(10) 
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^7 + 0(A a ), 



This means that A'" contributes mainly to Rayleigh scat- 
tering. Nevertheless, the term proportional to A a is im- 
portant for one-phonon Raman scattering. To first order 
in A„ we find that 



J 



( 3 "'l 



-a< 3 art° e A a ^ a 



gO a ,gl a {E e . g - AlflUJa - ftjJi + iT' e ) (E e>g + (1 - A 2 a ) fej B - fcj, + iT' e ) ' 

On the other hand, substituting (||) and (|io| ) in (^|) we obtain that, up to first order in the electron-phonon 
interaction, the B term is non-null only for one-phonon Raman processes. For Stokes processes at low temperature 
(i a = 0, j a = 1) we obtain 



(*"') 



/ go a ,gi a fr* (E e , g + (1 - A 2 a ) hw a - fkJi + iT' e ) (E e - E.) (E e>g - A 2 a hw a - huji + iTQ {E e - E s 



In the second term, one can exchange the indexes e and s. Next, under the condition \E e — E s \ ^> \Hu> a + i(T' e — T' s )\ 
one obtains 

g„A = -Wll g Wll s K e x a0 /V2 

)gO a ,gl a (£ e>9 + (1 - A*) hw a - fiw, + »T' e ) (E a , g - A 2 a ftj a - ft* + iT' s ) " 1 ' 



Noting that A 2 a <C 1, A a ftw a = — (e \He-l{o-)\ e), and we see that A'" + B'" for one-phonon emission reduces 
K e x a0 /V2 = (e \H E -l{o) \ s) (see Appendixes [A| and||), 
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TABLE II: Parameters used in the calculations When not 
indicated, the source is Ref. g for CdS and Ref. [| for CdSe. 

Parameter CdS CdSe CdSe (MBEMA) 



E g (eV) 

m e /m 

m h /rn 

7i 
72 

2m P 2 (eV) 

K 

V e (eV) 
V h (eV) 
u l (cm -1 ) 
ujt (cm -1 ) 
eo 



0L 



(10 -6 ) 
(meV) 



2.6 

0.18 

0.51 



21 c 

7.8 

2.5 

1.9 

305 

238 

8.7 d 

5.3 

2.68 

5 



1.865 

0.12 

0.45 



20 c 

9.53 

00 

00 

213 

165 

9.53 

5.72" 

1.58 

5 



1.841" 
0.13 6 

1.66 b 

0Al b 

20 c 

9.53 

0.6 b 

oo b 

213 

165 

9.53 

5.72 d 

1.58 

5 



a Ref. H 
b Ref. 24 
c Ref. \Z5 

^Calculated from the Lydanne-Sachs-Teller relation. 



to the same result that the TDPT. 



III. CALCULATION OF A'" AND B'" 

A. Effective Mass Approximation 

We use the exciton wave functions and the electron- 
phonon operator of Ref. ||(] to estimate the A and B terms 
in nanocrystals of several semiconductors. Using the Ta- 
ble I and the Effective Mass Approximation we have 



WT = 



where f 0e is the envelope overlap integral (o e being the 
set of envelope quantum numbers) 



\E> 0e (r,r)* d 3 r. 



p Sz j z -e is the bulk momentum matrix clement between 
the couple of bands to which the exciton e belongs. For 
the valence band J z = ±3/2, ±1/2 and for the conduction 
band s z = ±1/2. 

We evaluate the relative importance of the terms A and 
B for semiconductor nanocrystals. Assuming only one 
vibrational mode and focusing on one-phonon creation 
processes at low temperature (i. e. i — and j = 1), the 
resonant term in (^]) is reduced to 



(*") 



(gi\ev) (ev\gj) 



E 



E eig + (v — A 2 ) Hlulo -hioi+ iP 



-■ (12) 



The first term between parentheses is a band factor and 
is the responsible of the angular pattern of the scattered 
intensity. The summation over s z ,J z , is performed to 
take into account the degeneracy of conduction and hole 
bands, givingpl 



E 



(2J+1) 



(m Pf 1, 



(13) 



where J = 3/2 (1/2) if the upper valence band have Tg 
(r7) symmetry and P = —i (S \p x \ X) /mo- 

To evaluate the B term we also need to in- 
clude the electron-phonon matrix elements xoh se = 
V2 (e \He-l\ s) (see Appendix |b|). With these consid- 
erations we obtain 



E 



Ep 

foefos {e\H E - L \s) V2 




.!•() 



{gi\Q\ev) (ev\gj) + (gi\ev) (ev\Q\gj) 
E e , g + (v — A 2 ) Hlolo -hui + iT' e 



(14) 



We evaluate {c\He-l \ s) as in Ref. ^fj. Due to the 
Frohlich interaction cannot cause intersubband transi- 
tions, p* j . s = p* j . e and the band factor in (|l4| ) is 
the same as that in (|l^). Intersubband transitions may 
occur via deformation potential interaction, but these 
are usually negligible in polar materials. It must be no- 
ticed that it is not possible to consider this mechanism 
within the Albrecht's theory, as the existence of degener- 
ate hole states connected by the electron-phonon interac- 
tion means a breakdown of the adiabatic approximation, 
which is reflected in null denominators E e _ s in Eq. (|l4|). 
Nevertheless, TDPT can deal with it without trouble. 
From this result we conclude that the A and B terms 
are scalars (with small tensor corrections for B) and the 
integration over nanocrystal orientations and light polar- 
ization has no effect on the ratio between the A and B 
terms. 



B. Multiband effective mass theory 



In a multiband formalism, the essential effect of band 
mixing can be captured using-.the spherical approxima- 
tion for the hole HamiltoniancZl 



Hh 



7i 
2TOn 



p2 _!(p(2). j(2 ))) +m 



V(r) being the confinement potential, P^ 2 ) and J^ 2 ) are 
spherical rank tensors built from linear and angular mo- 
mentum operators, /1 = 272/71, and 72 and 71 are Lut- 
tinger parameters. 

Electron-hole pair states with well defined total 
(Bloch+orbital) angular momentum quantum numbers 
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M and M 7 can be obtained as 



l e ) — E ^nNlLsJfFMM 
n,N,l,L,f,F 



E ^ 



nNlLsJ fFMM z ) 
• > . Y / L.-.l fFMMt ifFfz^z I MM Z ) 



n,N,l,LJ,FJ z ,F z 



x |n/«//,) \NLJFF Z ) , 

where (fFf z F z \MM z ) is a Clebsch-Gordan coefficient. 
Lowercase (uppercase) letters denote electron (hole) 
quantum numbers, \nlsf f z ) and \NLJFF Z ) are electron 
and hole states with well defined total angular momen- 
tum, their wave functions given by 

(r\nlsff z ) = (lsl z s z \ff z )Rni(r)Y lh (0,<P) (r\aa s ) , 

lz,S z 



hole Bloch functions are related with the Tg (J = 3/2) 
electronic Bloch functions \J,J Z ) by the rule \JJ Z 



J, — J z ) (derived from the time-reversal opera- 
tion). Our | J J,) are |3/2,±3/2) = T(i/\/2)(X±iY) |±), 
and |3/2, ±1/2) = (t/V§) [2Z |±) =F (Jf ± iT) |t)]- 

Within this basis, in the strong confinement regime, 
the Coulomb interaction can be treated by direct diago- 
nalization of the Hamiltonian or even by simple pertvaj- 
bation theory. Using the theory of angular momentumEii, 
compact expressions for the matrix elements can be ob- 
tained. The dipole matrix elements gives 



2Peh 



Vi(-l) /+5/2 / 



(2/ + l)(2F+l) 



Ee,g V 3 

< ) 3 f 1, p\{S l , L + 8 l>L+2 ) J R^i L \r)R nl {r)r 2 dr, 



and 

(r\NLJFF z ) = E E ( KJL Mf.fz) 

K=L,L+2L Z> J Z 

xR$£\r)Y KL ,(0,<p) (r\JJ z ). 

In the above expression i? n ;(r) are the radial wave func- 
tions of a particle in a spherical box and RpjK^ (r) ar.c the 
solutions of the MBEMA equations given elsewhereEaEj, 
Yu z (9,ip) are the spherical harmonics^, (r|ss z ) are T 6 
Bloch function and (rl JJ Z ) are hole Bloch functions. The 



where R n i and R^^ are electron and hole radial func- 
tions, respectively, and 



e = k, e±i = =F 



i±jj 

V2 ' 



i, j and k being the unit vectors along the X-, Y- and 
Z-axis, respectively. 

In the A term we must sum over the degenerate e states 
with different M z , which turns out in a term proportional 
to the diagonal tensor 



y-yj 



V 4PW (2/ e + l)(2F e + l) f 1 3/2 1/2 

1 (>M,l((ll I ,l t + »I I ,L t +2) 7, 



E 2 

e,9 



R^y L '\r)R nJs (r)r 2 dr 



E 



(gi\ev){ev\gi) 



The B term for l p = phonons is given by 



S e . g + (v - A 2 ) Tkjlo ^ + iT' e 



/e+/a+3 



4P 2 e 2 ft 2 f 1 3/2 1/2 1 f 1 3/2 1/2 1 
ip =o ~ ^^S^^U /e F e /\Z S /. F s ft- 1 ' 

x V(2/ e + l)(2/ s + 1)(2F C + 1)(2F S + 1) (<^, ie + 5 lM2 ) {5 lstU + <5 is , is+2 ) 



x / R [ Nj;\r)RneiSr)r 2 dr J R^ (r)R nala (r)r 2 dr (i 

x V — ^\Q\ev) (ev\gj) + jgi\ev) jev\Q\gj) 
^ x (E e>s ) (E et3 + (v-i- A 2 )hoj LO -hu l + iT' e ) . 



H 



(i P =o) 

E-L 



IV. DISCUSSION 



r 



lated for typical nanocrystals 20 A in radius. The photon 
energies, corresponding to incoming resonance with the 
lower exciton level are also indicated in the table. The 



Table shows the Huang-Rhys parameter A 2 , the incomplete exciton confinement in CdS nanocrystals has 
A term absolute value, and the ratio \A"'/B"'\, calcu- 
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TABLE III: Numerical results for typical QD's 20 A in radius. 
The energies huii correspond to incoming resonance with the 
lower Raman active exciton level in each nanocrystal. 



Nanocrystal 


A 2 


huii (eV) 


\A"'/B"'\ 


\A'"\ 


(A 3 ) 


CdS 


0.08 


2.870 


13 


8.1 


X 


10 4 


CdS (R ef ) 


0.0013 


2.878 


8.4 


1.2 


X 


10 4 


CdSe 


0.0008 


2.592 


27 


7.9 


X 


10 3 


CdSe (MBEMA) 


0.2 


2.280 


32 


6.4 


X 


10 4 



been considered with two different models: (1) Finite 
band offsets V e and V^Ej, and (2) An effective radiusEl 
Moreover, for CdSe nanocrystals-we have also considered 
two models: (1) Effective radiusEl and (2) Multiband ef- 
fective mass approximation. (MBEMA) along with finite 
conduction band offset Y£p 

As it has been noticedEJEll, Raman scattering is quasi- 
forbidden in quantum dots in the strong confinement 
regime. The scattering is possible through a decompensa- 
tion between the electron and hole wave functions, which 
may appear due to (1) the Coulomb electron-hole interac- 
tion and difference between electron and hole masses, (2) 
difference in electron and hole confinement, (3) hole band 
mixing, (4) defects, and (5) non-adiabatic effectal-l The 
electron-hole decompensation also has a direct relation 
with the Huang-Rhys parameter. This one is larger in 
the case of CdS, where incomplete confinement have been 
considered, and in CdSe when band mixing and incom- 
plete electron confinement have been included. These 
ones are the cases where larger Raman polarizabilities 
are obtained. Notice that, due to the finite electron con- 
finement assumed, our Huang-Rhys parameter for CdSe 
nanocrystals within MBEMA is larger than other re- 
ported theoretical valuealj'EJ and is witix-the order of 
magnitude of the experimental valuesaornr l. 

In all the cases examined, the A term determines the 
Raman polarizability. This means that the interpreta- 
tion of one-phonon Raman cross section considering only 
the A term is consistent with the microscopic calculations 
using TDPT reported here and in Refs. How- 
ever, it is incorrect to fit the Huang-Rhys factor from the 
overtone to fundamental intensity ratios in multiphonon 
Raman spectra, as different scattering | chan . n els give sub- 
stantial contributions to the overtonestStSllH. 

Let us consider the participation of optical phonons 
with l p > in Raman scattering. These phonons connect 
degenerate band-mixed exciton states and breaks the 
adiabatic approximation. Hence, the Albrecht's theory 
cannot describe Raman scattering from these phonons. 
On the other hand, the TDPT, may deal with degen- 
erate states and non-adiabatic processes. Calculations 
in Ref. U indicate that the role of l p > phonons in 
one-phonon Raman spectra is to cause a small shoulder 
near the interface phonon frequency, being unimportant 
for the principal peak. Nevertheless, the electron-lattice 
interaction breaks the degeneracy of the exciton states, 
causing a redistribution of the exciton-phonon energy 



levels and possibly originate exciton-phonon complexes. 
Both these factors could substantially alter the predicted 
Raman spectra. This effect can be considered by TDPT 
of higher order in the electron-lattice interaction, or by a 
non-perturbative calculation of the exciton-phonon com- 
plexes. Research on this direction is presently in progress. 



V. CONCLUSIONS 

We have established the connection between the two 
theories more used for resonance Raman scattering in 
semiconductor nanocrystals: the Albrecht's theory based 
on the adiabatic approximation and the time dependent 
perturbative approach. In particular, we have shown 
that both theories are equivalent in the limit of weak 
electron-lattice interaction and when the resonant exci- 
ton level is non-degenerate or is a set of uncoupled degen- 
erate states. We have evaluated the relative importance 
of the Albrecht's A and B terms (the last one not dis- 
cussed in the literature) for CdSe and CdS nanocrystals, 
using different models for the electronic excitations. We 
have found that the A term is the leading coefficient in all 
the cases considered. Additionally, we have given the ex- 
pressions of the matrix elements of the electron-radiation 
and the electron-lattice interactions for a model of exci- 
ton considering the fourfold degeneracy of the Ts valence 
band. In the framework of this model, we have obtained 
a theoretical Huang-Rhys parameter within the order of 
magnitude of experimental values. 
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APPENDIX A: THE OSCILLATOR OFFSET 

The adiabatic vibrational eigenstates for the electronic 
state e are obtained from the Hamiltonian 

Hvibr = ( e \H-L + He-l \ e) 

= ]T {hw a b{b a + (e \H E _ L (a)\ e) (b a + b{ ) } . 

a 

We have considered the case of Hermitian He-l{o)- 
For non-Hermitian He-l{o) see Appendix |b|. The 
unitary transformation c a = b a — a a , with a a — 
— (e \He-l{<i) \ e) /hu> a diagonalizes the Hamiltonian 

H tibr = ^(clca - a 2 a ). 
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The new and old phonon coordinate operators are related 
by 



(c a + 4) = Q° - V2J —a„ 



2uj a - V w Q 

This relation identifies a a with the oscillator offset A a . 



in a semiconductor nanocrysta]0. Making the transfor- 
mation to coordinate and momentum operators 



the interaction operator becomes 



APPENDIX B: THE CONNECTION BETWEEN 
THE CLASSICAL AND THE QUANTUM 
ELECTRON-PHONON INTERACTION 

The Hamiltonian operator determining the electronic 
eigenstates can be expanded in Taylor series of the vibra- 
tional normal coordinates 

The second term of the above formula is the electron- 
phonon interaction. 



H E -l - 



OQa 



(Bl) 



The ouantum Frohlich-type interaction operator have the 
forrrO 

H e-l = E H E-L(o)k + H+_ L (a)bl 



where b a are annihilation operator of phonons in the 
normal modes a. H E _ L (a) are operators that act 
on the electronic coordinates, e.g. H EL (n,l,m) — 
$ n! ;(r)y; iirl (0, ip) for the one-electron-phonon interaction 



He-l = J2j^(H E _ L (a) + H+_ L (a))Q a 



H E _ L (a) - H+_ L {a)) P a . (B2) 



If H~_ L (a) = H+_ L (a) = H E - L {a) then (|Bl|) and (jB2|) 
are equivalent and 



h a 

I t„ <-■ 



dQ a 



2uj, 



^{s\H E . L (a)\e) 



When the vibrational modes in a nanostructure are de- 
scribed by complex fields there are modes for which 
H~ E _ L (a) ^ H E "_ L (a). In this case, thanks to time re- 
versal symmetry, complex modes are double degenerate 
and real fields can be obtained from the real and imag- 
inary parts of the complex fields, which correspond to 
real normal coordinates. The new matrix elements can 
be obtained from the complex matrix elements as 



H ( E l L (a) 



( S \H E _ L (a)+H+_ L (a)\e l 

. ( S \H E _ L (a)-H+_ L (a)\ 
V2 



(B3a) 



B3b) 



APPENDIX C: EXCITON PHONON MATRIX ELEMENTS WITH DEGENERATE BANDS 



The matrix elements of the interaction of band-mixed excitons with the optical phonons in a spherical nanocrystal 
can be calculated following the procedure outlined in Rcf. [l9| We obtained the expression 

(nN'l'L'sJf'F'M'M' z \H E _ L (n p ,l p ,m p )\ nNlLsJ fFMM z ) 
= (nNlLsJ fFMM z \H^_ L {n p , l p , m p )\ n N'l' L' sJ f F' M' M' z ) 

= ^ M '- K % ( ~U X M Z ) V(2M + l)(2M' + l)(2^ + l)/4 7r 
x {-S F , F >5 N , N ,5 LtL ,(-l) F + M+ f+f'+ s ^(2f> + 1)(2/ + 1)(2Z' + 1)(2/ + 1) 

x { f i v l s } { M f f f} ( 00)/ ^'W $ "^W^ r > 2dr 

+SfJ>5n,n>S l ,v{-l) M ' +f+J+2F Y, E y/( 2F " + l ^ 2F + + l ^ 2K + 1) 

K=L,L+2 K'=L',L'+2 

* { K K> ''j } { "f ")}{ 'S ) /^^H^IrA} ■ 
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where § np .i p {r) is the radial part of the optical phonons electrostatic potentialEl. The optical modes for l p > are 
described by complex fields, real field matrix elements can be obtained from Eq. (|B3|). 
For the Coulomb electron-hole interaction we obtained the expression 



n'N'l'L'sJf'F'M'M' z 



nNlLsJ fFMM z j = 5 M , M ,8 Mz , M ,{-l) s+J+2 f +F+F '+ M 

vW+W+ip+W+i) E { p 8 {> f i } { 

K,K',p ^ ' ^ 

x^(2l + l)(2V + l){2K + l)(2K' + l)^l P Q J) 



p / /' 1 f p f f> \ r p / /' 

J K' K \ \ M F' F 



V 


p / \ 




p 







i 










Rn'V {r e )Rni(r e )R%£> (r h )R y ^' (r h ) 



? ( F > L ) 



p+2 



dr„drh 



In the above expression, due to the properties of the 3j-symbols, the summation onp runs from max(|^ — l'\ , \K — K'\) 
to mm(l + l',K + K'). ^ 



* Electronic address: sariel@fT.oc.uh.cu 
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